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Abstract 

Gravitational perturbations in a Kerr black hole background can not be decomposed 
C"| ■ into simple tensor harmonics in the time domain. Here, we make the mode decom- 

position only in the azimuthal direction and discuss the resulting (2+l)-dimensional 
Klein-Gordon differential equation for scalar perturbations with a two dimensional 
t***'* ' Dirac's 5-function as a source representing a point particle orbiting a much larger 

black hole. To make this equation amenable for numerical integrations we explicitly 
remove analytically the singular behavior of the source and compute a global effective 
' source for the corresponding waveform. 

cr 

; 1 Introduction 

One of the main source targets of LISA is the gravitational waves generated by the inspiral of compact 
objects into massive black holes. For these Extreme Mass Ratio Inspirals (EMRI) we use the black 
hole perturbation approach to compute waveforms, where the compact object is approximated by a 
point particle orbiting a massive Kerr black hole. In order to obtain the precise theoretical gravitational 
waveform, we need to solve the self- force problem and problems in the second order perturbations pQ. 
■ In this paper, we focus on one aspect of the self-force problem, specifically to derive the retarded field 

of a point source. As a first step, we consider the Klein-Gordon equation in the Schwarzschild spacetime, 
but do not decompose it into spherical harmonics, in order to model perturbations in the more generic 
Kerr background. Recently, Barack and Golbourn [2] have discussed this equation in (2+l)-dimensions 
as derived by the mode decomposition in the azimuthal direction. Another treatment is proposed here 
to deal with this problem globally. We also note that there is a method to use a narrow Gaussian rather 
than a Dirac <5-function [3]. It is, however, difficult to ascertain the error introduced by smearing the 
^ ' particle and if this is accurate enough for self force computations. 

a: 

2 Formulation 

When we calculate the (2 + l)-dimensional (D) equation derived from the AD Klein-Gordon equation 
by the azimuthal mode decomposition, the resulting equation is not exactly same as the (2 + 1)D wave 
equation. By transforming the scalar field, we can obtain an equation which includes the (2 + 1)D 
dAlambertian and a remainder. Then, we remove the 2D (^-function in the source term. 
In order to do so, we consider the Schwarzschild metric in the isotropic coordinates, 

ds 2 = -(2p - M) 2 / {2p + M) 2 dt 2 + (1 + M/{2p) f (dp 2 + p 2 (d9 2 + sin 2 6d<p 2 )) . (1) 
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This radial coordinate is related to that of the Schwarzschild, r, as p = (r — M + \/r 2 — 2Mr) /2. In the 
above coordinates, the Klein-Gordon equation with a point source becomes 

— -odt + — j d l + + £ — t d 2 + cot 6d g + — ^-i 

(2p-M) 2 (2p+M) 4 p (2p~M)(2p + M)° (2p + M) 4 \ sin 2 6 

= r dT 64A(t-t,(rM P -p (r)M( fl - g ,(r) W -0,(r)) } 
VJ V-oo (2p-M)(2p + M)5 s in0 w 
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Here, we use the azimuthal mode decomposition, tl>(t,r,9,(p) — X)m=-oo VVn^ P, 9)e lm ^ , and transform 
the field tp m as tp m = 2yJ p/[(2 p + M) (2p — M) sin#] x-m- Then, Xm satisfies the following equation. 



16(2p-M)V P P P 2 V sin 2 0V"" 4 (2p + M) 2 (2p — M) 2 

i-co v /(2p-M)(2p + M)psin0 W 

Next, we use a new time coordinate which is defined by T = f* dt 4 (2p z (t) — M)p z (t) 2 / '(2p z (t) + M) 3 , 
where p z is obtained by solving a geodesic equation. Using this, we derive an equation which can be 
divided into the (2 + 1)1? d'Alambertian of the flat case n( 2+1 ) and a remainder. 

C m Xm (T, p, 9) = (n( 2+1 ) + CZ m ) Xm(T, p, 9) = S rn (T, p, 9) ; 

□ (2+1) = -a 2 +a 2 + (i/p)a p + (i/p 2 )a 2 , 

1 _ {2 Pz (T)-M) 2 p z {T)\2p + Mf \ 2 

(2 Pz (T) + M)6(2p-M)V ) T 
2{Ap z (T) - M)(2 Pz (T) - M){2p + Mfp z (TfM f dp z (T) 



(2p z (T) + My(2p~ M) 2 p 4 V dT 

2 1 (V + M 2 ) 2 - 16pMcos 2 9" 



p 2 sm 2 9\ 4 (2p + M) 2 (2p-M) 2 

S m (T,p,9) = 2 r dT ^m^!kzmm^^e-^ (4) 
J-oo y/(2p-M)(2p + M)psmO 

To remove the (^-function in the source term, we set 

Xm (T,r,9) = X S m (T,r,9)+x^ m (T,r,9), (5) 

where we define the new functions, Xm an d X™ as calculated from 

□ (2+1) X*(2>,0) = S m (T,r,9); L mX ™ n {T^9) = ^ m X S m (T,r,9) = S^(T,p,9). (6) 

This decomposition of Xm does not have any physical- meaning, i.e., Xm 1S n °t identified as the singular 
part to be removed in the self-force calculation. Note that £„ m includes a second-order derivative. But, 
since the factor of d\ is zero at the particle location, the singular behavior of the effective source S^xi^ 
weakens. The derivation of the singular field \fn can ^ e performed through the Green function, 

G(T,x;T',x') = 9((T - T') - |x - ^\)/{2^(T - T' f - |x - x'| 2 ) , (7) 
where |x — x'| = (p 2 + p' 2 — 2 p p' cos(9 — 9 1 )) 1 / 2 and Xm IS calculated by the following integral 

xi(T,p,9) = JdT'p'dp'd9'G(T,x;T', X ')S m (T',p',6'). (8) 

3 Circular Orbit Case 



We consider a particle in circular orbit, z a (r) = {u*t, ro, n/2, u^t}, where w* = \/ r o/( r o — 3M) and 
%& = ^/M/[ro(ro — 3M)]. Here, we focus on the m modes; the to = mode can be dealt with by 
the same treatment. 

3.1 Singular field 

First, the relationship between the new time coordinate T and the Schwarzschild time t is obtained 
analytically as T = 4 (2p - M)pg/(2p + M) 3 t where p = l/2(r - M + y^o = 2Mr ). Note that in 
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general, for non circular orbits, we need a numerical integration to derive this relationship. From Eq. ([5]), 
the singular field is derived as 

X S m (t,r,e) = 1 - 2 ^~° H^( A { jj P0 + ^ 2 mnix-x.lV-^, (9) 

where f2 = u^/u* and the spatial difference is given by |x — x z | = ^/ p 2 + p 2 — 2 p po sin and -ffg 1 ^ is the 
Hankel function. The local behavior of the above solution near the particle location is obtained as 

Xl(t, r, 9) ~ X i L (t, r, 9) = -±- - = 2 ^ = ln f n»fl|x - x z |) . (10) 



2tt u* A /(2p + M)(2A)- M) V 4 (2Po-M)Po 

Therefore, 5™ (i, p, 9) shown by the dashed green curve in Fig.[TJ is singular at the particle location. In 
order to perform the numerical integration with higher accuracy, it is convenient to regularize the source 
term to be at least C° at the particle location. 



3.2 Local behavior 

When we write the singular field as Xm = Xm L + Xmi Xm 1S finite at the particle location. Then, the 
effective source in Eq. ([6]) becomes 

S^Ht, P ,9) = * ,. 2 ^ lnr . (2 /° +M)3 . m01x-x z lV-^ 



2tt u t y/(2p Q + M){2p - M) \A{2p -M)p 2 
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1 (4p 2 + M 2 ) 2 - 16pM cos 2 9 



\ p 2 sm 2 9\ 4 (2p + M) 2 (2p- M) 2 

1 ' \ (2p - Af) 2 p4 (2p - M) 2 p± ) J Lm Xm[h ' ' ' 

Note that the third line of the R.H.S. is at least C° at the location of the particle. 
To remove the logarithmic divergence in the source, we introduce 



X^ S (t,P,0) = |x - x z | 2 In ( 4 [ 2 /; Q + ^ m O,x - x z |) p ^/ 2 (2 p - Mf e* 

64 m 2 p 4 - 32 m 2 p 2 M 2 + 4 m 2 M 4 + 16 cos 2 0p 2 Af 2 - (4p 2 + M 2 ) 2 



(12) 



ttu* (2 p + M) 5/2 (2 p - A/) 11/2 p 11 sin 2 9 

Using this regularization function, we obtain a source S^f 9 ' 7 f° r the function £™ — x™' S as 

^• 7 (t,p,0) - S^f\t,p,9)-C mX r ^ m ' S (t,p,9). (13) 

The local behavior of S 1 ^ 9 ' 1 , which is shown by the solid green curve in Fig. [1] (b), is an "ccln|x| for 
x — ► 0" type, i.e., C° around the particle location. 



3.3 Boundary behaviors 

We now focus on the behaviors of the source term at the two boundaries, i.e at the horizon of the large 
hole and spatial infinity. The source for a final regularized function Xm 9 must go like 0(p~ 2 ) for p — > oo 
in the case of the m = mode and 0(p~ 3 / 2 ) for the m ^ mode because of integrability conditions. More 
precisely, the source for the regularized function of ip m derived by numerical calculations has a factor 
- l/p 1/2 . For p -► M/2, the source should be zero, i.e., the behavior should be a power of (p — M/2) 
greater than 1/2. To regularize the source at the boundaries, we note that the source contribution from 
Xm m ' S 1S wen behaved. This means that the ill behaviors of the source arise from x m - Therefore, it is 
convenient to use asymptotic behaviors of X fn (and some correction factor) for regularization. 
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First, for the regularization near the horizon, we use the regularization function x m which is too long 
to be shown here. Then, the source for the function Xm™ ~ X™' S — Xm becomes 

S% 9 > h (t,p,6) = 5^ 7 (t,p,0)-£ mX ^(<,p,0), (14) 

This Sm 9 ' h is shown by the solid black curve in Fig. [Hand behaves as (9(p -1 / 2 ) for large p. To regularize 
it, we use the regularization function, 




The final source for the regularized function Xm 9 = X™ — X™' S — X m ~ Xm becomes 

s%>>f(t,p,o) = s^ h {t, P ,e)-jC mX ^(t,P,e), (ie) 

which is used in the numerical calculation. This S 1 ^ 9 is shown by the dashed black curve in Fig. [T] and 
behaves like 0(p~ 3 / 2 ) x (an oscillation factor with respect to p) for large p. 
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Figure 1: Plot for the m = 1 mode of S m with respect to p. S[ e ^ , S[ e9,1 , S'[ e9 '' i and S[ e9 ^ are shown 
by the dashed green, solid green, solid black and dashed black curve, respectively. 



4 Discussion 

In this paper, we obtained the regularized effective source which is C° at the location of the particle, and 
0(p — M/2) near the horizon. The behavior at infinity is 0(p~ 3 ^ 2 ) x (an oscillation factor with respect 
to p) which allows straightforward numerical integration. 

When we consider the extension of this formulation to the Kerr background case, we can also extract 
a similar differential operator to that of Eq. (j4|). In the case of gravitational perturbations, we have ten 
field equations for the linear perturbation in the Lorenz gauge. (See [4].) The same treatment discussed 
in this paper is applicable to those equations. 
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